Lane-formation vs. cluster-formation in two dimensional square-shoulder systems: 

A genetic algorithm approach 
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Introducing genetic algorithms as a reliable and efficient tool to find ordered equilibrium struc- 
tures, we predict minimum energy configurations of the square shoulder system for different values 
of corona width A. Varying systematically the pressure for different values of A we obtain complete 
sequences of minimum energy configurations which provide a deeper understanding of the system's 
strategies to arrange particles in an energetically optimized fashion, leading to the competing self- 
assembly scenarios of cluster-formation vs. lane-formation. 

PACS numbers: 



The ability of colloidal dispersions to self-organize in a 
surprisingly large variety of ordered structures markedly 
distinguishes soft matter from hard (i.e., atomic) ma- 
terials. The spectrum of ordered particle arrangements 
encountered in soft matter systems does not only con- 
tain low-symmetry, non-close-packed structures, such as 
bco, diamond, or the A15 lattices [H, 0, 0, 0, [H; s °ft 
matter particles are able to self-organize in considerably 
more complex ways, forming thereby micellar and inverse 
micellar structures @, 0], cluster phases 01, chain-like 
or layered arrangements @, 0, EH EH EH EH , or gyroid 
phases EE EH > J ust to name a few of them. 

In the investigation of self-organizing phenomena in 
soft matter systems, theoreticians have developed reliable 
tools. On one side efficient and accurate coarse-graining 
procedures have been developed which, by averaging over 
the large number of degrees of freedom of the solvent 
particles |l7|, [l8| , lead to effective potentials between the 
colloidal particles. On the other side there are statistical 
mechanics based concepts which provide on the basis of 
these interactions reliable information on the thermody- 
namic properties that finally lead to the phase diagram 
of the system at hand. Among those are density func- 
tional theor y Il9| . liquid state theories |2p| . or computer 
simulations [21] . 

What is badly missing among the theoreticians' tools 
is a reliable way of how to predict the ordered equilib- 
rium structures for a given system. In hard materials 
one can safely rely on experience, intuition, or plausible 
arguments: a set of pre-selected candidate structures is 
chosen, comprising the usual suspects (such as fee, bec, 
etc.) and at a given state point the lattice with the lowest 
(free) energy is the stable one. In soft matter, however, 
the situation is entirely different: as a consequence of the 
rich wealth of emerging equilibrium structures a conven- 
tional approach that is based on a biased pre-selection 
process is bound to fail. Computer simulations are not 
very helpful cither, since they are time-consuming and 
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risk to be trapped in local energetic minima due to the 
rough and complex energy surface. 

In this contribution we present an alternative strategy 
which we believe to be very powerful in the search of equi- 
librium structures and which does not share the deficien- 
cies of the conventional approaches. We use genetic al- 
gorithms (GAs), i.e., optimization strategies, that adopt 
features of evolutionary processes as key elements to find 
the optimal solution for a problem [12] • In contrast to the 
conventional methods, GAs allow for searching basically 
among all lattice structures in a parameter-free and unbi- 
ased way. Despite the fact that GAs have been proposed 
several decades ago [HJ, their usefulness in problems of 
condensed matter theor y h as been acknowledged only in 
recent years (see, e.g., [23l . H3])- In the present contri- 
bution we demonstrate the power of this approach for a 
two-dimensional (2D) system where the particles inter- 
act via spherically symmetric square-shoulder potentials. 
With the help of GAs we identify minimum energy config- 
urations (MECs), or, equivalently the zero-temperature 
phase diagram. 

GAs mimic certain principles and processes known 
from natural evolution like mutation, mating and 'sur- 
vival of the fittest' to find the global extremum in the 
function to optimise. Because of their particular search 
strategies they are able to investigate large areas in search 
space and to concentrate at the same time their efforts 
on the most promising regions. It is in particular this 
global scope which has made GAs a highly appreciated 
tool in many fields. 

The basic unit of a GA is a so-called individual, X, 
which represents one possible solution of the problem, 
encoded in a string of genes. Each gene can take on 
values out of a chosen alphabet. The encoding of the 
candidate to the individual is a very delicate task that 
has major influence on the performance of the GA. A set 
of individuals is called generation. 

With these definitions we can sketch the work-flow 
of a GA: first a starting generation is created at ran- 
dom, whose individuals consist of arbitrary sequences of 
genes. Each of the individuals T is assigned a 'fitness- 
value' through a problem-specific function f(T), assess- 
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ing the quality of the solution represented by the indi- 
vidual. A higher fitness-value marks a better solution in 
this evaluation process. According to their fitness, par- 
ent individuals are selected for reproduction to generate 
new individuals by simple operations like recombination 
and mutation, for which many different methods have 
been proposed in literature. This circle, i.e., evaluation- 
selection-recombination-mutation, is iterated until a suf- 
ficiently large number of generations has been created. 
The individual with the absolutely highest fitness value 
is taken as the final solution. Since GAs do not converge 
exactly to the global minimum, additional refinement of 
the proposed solution is required. 

In an effort to adapt this general algorithm to our par- 
ticular problem, we first have to find a suitable encoding 
of the candidate structures to the individuals. To this 
end we have utilized two different strategies: (i) a simple 
lattice paramctrisation where the lattice vectors and po- 
sitions of eventual basis particles are translated via the 
binary alphabet to strings of Os and Is; (ii) a more refined 
strategy, which is specially adapted to systems for which 
cluster phases are to be expected, has been developed: 
in this 'cluster-biased' version the lattice sites are popu- 
lated by ordered (i.e. regular) dimers or trimers. Hence 
individuals include additional information describing the 
orientation and the distance of particles in the cluster. In 
this way more particles per unit cell can be considered, re- 
ducing at the same time the number of parameters. This 
leads to a considerably enhanced performance of the GA 
for complex lattices. In addition, for hard-core particles, 
particular care is in order to prevent the creation of in- 
dividuals that represent unphysical configurations with 
overlapping cores. 

Working at fixed particle number N, pressure P, and 
temperature T, we have decided for the standard form of 
/(J), i.e., 

/(T)=exp{-[G(X)-G ]/G } . (1) 

At T = 0, the Gibbs free energy, G, reduces to G — 
U+PV, V being the volume. Go stands for the Gibbs free 
energy of a reference structure. For further conceptual 
and numerical details concerning the GA we refer to [25|, 

m. 

The interaction potential of the square-shoulder sys- 
tem consists of an impenetrable core of diameter a with 
an adjacent step-shaped, repulsive corona (with range 
Act), i.e., 

$(»") = \ \ , (2) 

w \ e CT<r<Aer' w 

e being the height of the shoulder. The search for MECs 
for this particular system represents undoubtedly the 
most stringent test for the GA: with its flat plateau and 
as sharp a cutoff as possible the MECs can be easily clas- 
sified by the number of overlapping coronas which makes 
the potential a quintessential [2j test system. The fact 
that we restrict ourselves to a 2D system does not rep- 
resent a limitation at all. On the contrary, the particle 



arrangements can be visualized in a very convenient way 
that makes it much easier to understand the sequence of 
MECs as the pressure is increased. Objections that the 
system is oversimplified and only of academic interest can 
easily be refuted: there exist a number of realistic soft 
systems with a core-corona-architecture that can be de- 
scribed using a square-shoulder interaction fl3j . Among 
these are, for instance, colloidal particles with block- 
copolymers grafted to their surface where self-consistent 
field calculations lead to effective interactions that closely 
resemble $(r) given in (2) (27| . 

The search for MECs for the square-shoulder poten- 
tial has been pioneered by Jagla [§, [l(| and was carried 
on in later work by Malescio and Pellicane [12, EH using 
Monte Carlo simulations and geometrical considerations. 
Further theoretical work on this system was presented in 
1\. The results are indeed remarkable: a large variety 
of MECs has been identified where - despite the radial 
symmetry of $(r) - particles often arrange in asymmetric 
structures, forming thereby lanes or ring-like structures. 
However, the authors of above mentioned studies raise 
doubts themselves that their set of MECs might not be 
complete. Introducing GAs as a novel technique to this 
particular problem, we will present in the following se- 
quences of MECs that show a considerably larger variety 
than the ones identified up to now and which can eas- 
ily be understood via energetic arguments. Although we 
cannot provide a rigorous proof there is evidence that 
these sequences are complete. 

Our results offer - depending on the range of the shoul- 
der A - a deeper insight that the square-shoulder system 
develops different strategies to form MECs. Among the 
systems we have investigated we present in the follow- 
ing results for three different values for A, correspond- 
ing to a small (A = 1.5), an intermediate (A — 4.5), 
and a large shoulder range (A = 10). Standard re- 
duced units are used: P* = Pa 2 /e, U* = U/(Ne) and 
G* = G/(Ne) = U* + P7(//ct 2 ), T] = N/V being the 
number density. 

We start with A = 1.5 and show in Fig. Q]the MECs 
proposed by the GA, while Fig. [5] displays the corre- 
sponding thermodynamic properties G* and U* as func- 
tions of P*. At very low pressure, particles populate an 
ideal hexagonal lattice, thus avoiding overlapping coro- 
nas. Upon compression, the system must pay in some 
form tribute to the reduced space in terms of an energy 
penalty, i.e. via a first overlap of shoulders. Obviously, 
for this A-value, the formation of lanes is energetically 
the best solution. Along these lanes, particles are in di- 
rect contact (and, consequently have overlapping coro- 
nas), forming a one-dimensional close-packed arrange- 
ment. Parallel lanes, however, try to avoid corona overlap 
and the shoulder width A serves as a spacer (see magnified 
view). As the pressure is further increased, new strate- 
gies are required. While particles still prefer alignment 
along lanes their internal arrangement is modified: rather 
than forming straight lines, the lanes are zig-zag shaped 
due to energetic reasons, which is a compromise between 
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FIG. 1: MECs for the square-shoulder system of shoulder 
range A = 1.5. Configurations correspond (from left to right 
and from top to bottom) to pressure values indicated in Fig. 
|2] by vertical arrows. 
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FIG. 2: G* (black line) and U* (grey line) as functions of 
P*/X 3 for the square-shoulder system with A = 1.5. Verti- 
cal arrows indicate MECs depicted in Fig. [1] Broken lines 
represent limiting cases of MECs (see text). 



the reduced available space and the imminent energetic 
penalty due to additional corona overlaps. Neighbouring 
lanes are arranged in such a way, that each particle is now 
in direct contact with three other ones. Alternatively, the 
staggered lanes can also be viewed as a ring-like struc- 
ture: six particles form elongated rings where A serves 
again as a spacer, fixing the width of the cage. In the end, 
further compression causes the system to collapse into the 
close-packed hexagonal structure where each particle is 
in direct contact with six neighbours. 

The functional form of G* and U* can nicely be under- 
stood via the following, general thermodynamic consid- 
erations (cf. Fig. [2]): For a given particle configuration, 
characterized by the number of overlapping coronas, U* 
is constant; thus G* = U* + P* /(ija 2 ) is a straight line 
as a function of P* with slope 1/rj. Two limiting cases 



(indicated as broken lines in Fig. [2] and characterized by 
slopes 1 / rj m m and 1 / ?7 m ax) can easily be identified: an ar- 
rangement where disks of diameter A form a close-packed 
structure (with ^ m i n and U* = 0) and the hexagonal lat- 
tice where the particles' hard cores form a close-packed 
structure (with ?7 max and U* = U* lax = 3). All other 
MECs identified by the GA are located on lines of slope 
I/77, with l/r? m ax < I/77 < l/?7niinj an d obviously have 
lower G*-values. Thus G* = G*(P*) becomes a sequence 
of intersecting straight lines, representing a considerable 
help in analysing the data. Cross-over points between 
two competing structures where GAs tend to have con- 
vergence problems [25[ can be determined exactly as the 
intersection of two straight lines. The energy levels U* 
(also depicted in Fig. [J) that characterize the MECs are 
of course rational numbers, i.e., number of overlapping 
coronas divided by the number of particles per unit cell. 
Finally, all calculations have been performed on rather 
fine pressure grids (AP* ~ 0.2) so we are confident that 
the presented sequences of MECs are complete. 

As we proceed to A = 4.5 the systems develops com- 
pletely different strategies to form MECs as the pressure 
is increased; they are depicted in Fig. [H while Fig. 2] dis- 
plays the thermodynamic properties. The hexagonal pat- 
tern imposed by the non-overlapping coronas (observed 
at extremely low pressure- values and not displayed here) 
is soon superseded by a novel strategy, namely cluster 
formation. At low pressure, dimers are formed which 
populate the sites of a distorted hexagonal lattice. Upon 
further compression, these aggregates become larger un- 
til they reach the size of six particles. The degree of 
distortion of the underlying hexagonal lattice is imposed 
by the shape of the clusters: therefore, the trimers, which 
have nearly circular shape, sit on a nearly perfect hexag- 
onal lattice, while for elongated hexamers the structure is 
strongly distorted. As the system is further compressed, 
formation of clusters is obviously energetically less at- 
tractive and lane-formation sets in. With increasing pres- 
sure the structure of the MECs becomes more complex. 
In the beginning each lane is built up by a linear se- 
quence of dimers and thus closely resembles the first lane- 
scenario observed for A = 1.5. At higher pressure- values, 
however, the system forms parallel lanes which are now 
built up by larger clusters. The increasing complexity of 
this inner structure makes simple energetic explanations 
in terms of overlapping coronas impossible. However, the 
grey shades in Fig. Ogive evidence that the formation of 
clustered lanes is an efficient strategy to avoid an overlap 
between neighbouring lanes, which represents a consider- 
ably higher energy penalty. The tendency to form paral- 
lel lanes is maintained until the system finally collapses 
into the close-packed hexagonal structure. The consid- 
erably richer wealth of MECs encountered for A = 4.5 is 
reflected by the large number of energy levels in the plot 
U* vs. P* and the number of intersecting straight line 
segments in the diagram of G* vs. P*. 

Finally, for A = 10, the strategies of the system to form 
MECs seem at first sight similar to the previous case: 
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FIG. 3: MECs for the square-shoulder system of shoulder 
range A = 4.5. Configurations correspond (from left to right 
and from top to bottom) to pressure values indicated in Fig. 
H]by vertical arrows. 
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FIG. 4: G* (black line) and U* (grey line) as functions of 
P*/\ 3 for the square-shoulder system with A = 4.5. Verti- 
cal arrows indicate MECs depicted in Fig. [3] Broken lines 
represent limiting cases of MECs (see text). 



again, the system prefers formation of clusters at low 
pressure, which are located on slightly distorted hexag- 
onal lattices. However, we observe that inside a cluster 
the cores of the particles are sometimes arranged in a dis- 
ordered fashion, while for A = 4.5 only clusters with an 
ordered internal particle arrangement occur. The strat- 
egy is obviously the following: once A is sufficiently large 
so that cluster formation is supported, the system tries 
to arrange particles so that the shape of the cluster be- 
comes as circular as possible. This in turn guarantees 
that the underlying structure is close to the energetically 
most favourable hexagonal lattice. For A = 4.5, where 
the core region still represents a considerable fraction of 
the particle diameter, the system has to proceed rather 
carefully to fulfill this requirement, leading to the or- 
dered arrangements of the cores. For A = 10, however, 
the core region is nearly negligible with respect to the 
core width. Now both regular and irregular particle ar- 
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FIG. 5: MECs for the square-shoulder system of shoulder 
range A = 10. Configurations correspond (from left to right 
and from top to bottom) to pressure values indicated in Fig. 
H3 by vertical arrows. 
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FIG. 6: G* (black line) and U* (grey line) as functions of 
P*/\ 3 for the square-shoulder system with A = 10. Verti- 
cal arrows indicate MECs depicted in Fig. [1] Broken lines 
represent limiting cases of MECs (see text). 



rangements inside the core can lead to circular-shaped 
clusters of the same size, having practically the same G- 
value. Again, at higher pressure values lanes with an 
increasingly complex inner structure are formed. Some 
of them (for example, the last one displayed in Fig. [5]) 
might be comparable to the Bernal spirals observed ex- 
perimentally in three dimensional colloidal systems [28j ] . 
Lane formation persists until the system finally collapses 
into the close-packed hexagonal structure. If the square- 
shoulder particles still represent a reasonable model for 
macromolecules at such high pressure values is question- 
able. 

Let us come back to the disordered clusters. If we ex- 
trapolate these results to even larger values of A - or, 
equivalently, to a vanishing core - and briefly switch to 
three dimensions, we arrive at another soft matter model 
system that has been studied in detail in literature: the 
penetrable sphere model (PSM) (29|. Among the remark- 
able, well-documented features of this system is its ability 
to form cluster phases: at sufficiently high densities clus- 
ters of overlapping particles populate the lattice sites of 
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a regular fee lattice. Detailed simulations of a closely re- 
lated system that also shows clustering Q have revealed 
that the internal structure of the clusters is completely 
random. These observations are consistent with our re- 
sults in the sense, that an increasing corona width (or, 
equivalently, a vanishing core) favours formation of dis- 
ordered clusters of particles, which in turn populate sites 
of regular lattices. Obviously first precursors of this phe- 
nomenon occur at A = 10 [3fJ. 

Our results provide an explanation for the system's 
strategy to form MECs: if A is small, formation of larger 
clusters is prohibited due to geometric reasons, hence 
particles arrange in lanes of variable shape or in con- 
nected structures. However, as soon as the range of the 
corona is sufficiently large, the formation of clusters sets 
in and they dominate the low-pressure regime. Depend- 
ing on the corona range, the particles group in ordered or 
disordered clusters, keeping the cluster's shape as disk- 
like as and the underlying structure as close to the ideal 
hexagonal lattice as possible. For systems that exhibit 
clustering, lane formation sets in at higher pressure val- 
ues. The lanes themselves are built up by clusters, so one 
is led to interpret the crossover to the lane regime as a 



structural change of the cluster crystal. 

The results of our investigations are, on one side, of 
rather basic relevance. The simple form of the poten- 
tial allows to understand the full sequence of MECs on 
an energetic level as the pressure is increased and thus 
provides deeper insight under which conditions the sys- 
tem prefers to form either lanes or clusters. But also 
from a more applied point of view, these results are of 
relevance. Having established GAs as a suitable and re- 
liable tool for finding ordered equilibrium structures, the 
strategies of self-assembly of systems with a substantially 
more complex effective interparticle interaction can be 
better understood. Such knowledge is certainly of rele- 
vance in technological applications, such as nanolithog- 
raphy or nanoelectricity. 
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